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ABSTRACT 

Turbulence models are developed by supplementing the renormalization group (RNG) 
approach of Yakhot k Orszag with scale expansions for the Reynolds stress and production 
of dissipation terms. The additional expansion parameter [rj = SK /£ ) is the ratio of the 
turbulent to mean strain time scale. While low-order expansions appear to provide an 
adequate description for the Reynolds stress, no finite truncation of the expansion for the 
production of dissipation term in powers of 77 suffices - terms of all orders must be retained. 
Based on these ideas, a new two-equation model and Reynolds stress transport model are 
developed for turbulent shear flows. The models are tested for homogeneous shear flow 
and flow over a backward facing step. Comparisons between the model predictions and 
experimental data are excellent. 
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I. Introduction 


According to the Kolmogorov theory 1 of turbulence, the dynamics of velocity fluctuations 
v(k,u>) at the scale / = 2ir/k depend only on the mean rate of dissipation of kinetic energy 
£ and the length scale l. Neither the integral (L) nor dissipation (Id = 2ir/kd) scales enter 
the dynamics of the inertial range of scales where Id <C l L. Based on this theory, the 
velocity of eddies of size l scales as £ ^ Z 1 ^ 3 so that the characteristic (turnover) time of these 
eddies is 

T(i)~-L^r U 3 p> 3 . (i) 

v(l) 

The velocity correlation function can then be represented in the scaling form 

< v<(k,u;)v,(k',ty) >~ k~ l3/3 <j> ^i/3 ~ fc2/3 ) - C ( k > u ) ( 2 ) 

where the scaling function (f>( x) is to be determined from other considerations. The energy 
spectrum E(k ) is of the form 

E(k) ~k 2 J C(k,uj)du = C K £ 2/3 k- 5/3 (3) 

where Ck is the so called Kolmogorov constant. These results also yield the effective (tur- 
bulent) viscosity vj at the scale 1 = 2i r/k: 

ux(k) ~ Vil ~ . (4) 

This effective viscosity plays a profound role in turbulence modelling. For example, 
it is quite reasonable to assume that the large-scale properties of the flow are governed 
by effective equations of motion similar to the Navier-Stokes equations with the effects of 
strong interactions between the velocity fluctuations taken into account through an effective 
viscosity (see Yakhot and Orszag 2 ) 

u T ~ r 1/3 L 4/3 . (5) 

However, the strict notion of eddy viscosity requires the existence of a small parameter 

l/L < 1 which is absent in turbulent flow. Still, the eddy viscosity concept proves to be 
extremely useful, working much better than expected. This situation is not unique: in a 
fluid close to thermal equilibrium, the molecular viscosity representation is very accurate 
even when the mean free path A is not that small (A/ L — 1). 

Many years before Kolmogorov’s work, Osborne Reynolds realized that turbulent flow was 
a random system which had to be treated using statistical methods. The work of Reynolds 
was similar in concept to the work of Boltzmann, Gibbs and others who formulated the 
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kinetic theory of gases. It was not an accident that by the early 20th century the ideas of 
kinetic theory were adapted to describe turbulence with turbulent eddies as the molecules 
or building blocks of this “gas.” By analogy with kinetic theory the turbulent viscosity is 
taken to be: 

V' T — ’U'rmsL* ( 6 ) 

Combining (6) with (5) and using £ ~ u^ m JL (so that, in the absence of turbulence pro- 
duction, the turbulence is damped in one large-eddy turnover time L/u rm s), we obtain the 
K — £ model: 

K 2 

VT = C/i-gT (7) 

where K — \uiul and C is a constant. The relation (7) is more convenient than (6) since it 
expresses the turbulent viscosity through the directly measurable quantities K and £. 

To implement (7) we need equations of motion for K, £ and the velocity field v: 


v • Vv - - Vp + i/ 0 V 2 v 
— + V • VI< = -£ - Vi(vip) + v 0 V 2 K 

A 

— + v • V£ = -Ti -T 2 — 2vo(V j v i )(V i V J p) + v 0 V 2 £ 
where p is the pressure and 

Ti = 2 vq (V J u,)(VjU;)(V i i;,-) 

T 2 = 2vq(V jV iVi) 2 . 


( 8 ) 

( 9 ) 

( 10 ) 


( 11 ) 


Equations (8)-(ll) must be solved subject to initial and boundary conditions as well as the 
incompressibility constraint V • v = 0. In order to derive equations for the mean values of K 
and £, equations (8)-( 11) must also be averaged over an ensemble of the fluctuating part of 
the velocity field v. In the absence of a systematic averaging procedure, turbulence modellers 
have used intuitive reasoning combined with experimental data, tensorial and dimensional 
invariance, as well as scaling arguments based on the Kolmogorov theory 3 ”" 6 . 

In the simple case of decaying turbulence with no mean velocity field, the equation of 
motion for the mean kinetic energy is usually written as: 


OK 

dt 




( 12 ) 


where v = vq + vj is the total viscosity and an is a proportionality constant. The derivation 
of the equation for £ is more difficult since it is easily shown that both T\ and T 2 in (11) are 
0(Re t ) where Re t = K /v q£ = 0(u rms L/ v 0 ) is the turbulence Reynolds number. However, 
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the following bold assumption (cf. Tennekes and Lumley 7 ) is fruitful: it is assumed that the 
singular contributions to T\ and T 2 cancel each other in the limit as Ret * 00 and the 
remaining 0(1) part can be written as 

2\ + T 2 = C £ S= (13) 

K 


where C £ 2 = 0(1). Thus, in decaying turbulence: 


d£ 

at 


£ 2 a ( os' 

= ■ Cc2 t + (r^ 


(14) 


where a £ = 0(1). The assumption (13) was not rigorously justified but it has led to simple 
equations useful for practical calculations. In the case of homogeneous isotropic turbulence 
the diffusion terms in (12) and (14) disappear yielding 


K oc r 7 


(15) 


where 7 = 1 f(C £2 - 1). The fact that the coefficient C £2 determines the power law of 
turbulence decay demonstrates the significant role of dimensionless constants in the theory 
of turbulence. Indeed, even a small error in C £2 is amplified in the limit as i —* 00 . The first 
term on the right side of (14) has a simple interpretation in terms of a relaxation time: 

t^__£_ 

T ~ T 


where — ^ 

T = J 

is the only turbulence time-scale that one can construct from the parameters of the problem. 

In flows with non-zero values of the mean strain 5,j = |(§~^ + 9^f)> velocity field in 
(8)-(ll) can be decomposed into mean (U) and fluctuating (u) components. In this case the 
equations of motion are typically taken to be of the form. 

dK — ^ d ( dl< 

- + U.VX = -r,5,-£ + -U^ 



dS 


s 2 


d 





where = uju] is the Reynolds stress tensor and the coefficients C £ 1 and C £2 are constants 
that are usually determined from benchmark experiments. The £ - production term in (16) 
comes from the contributions to the £ - equation (10) of the type: 


dui duj 
2Vo dxt dx, ij 
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which are modeled as 


CdTijSij/T 


in the relaxation time approximation. 

There are also dynamical constraints that a consistent turbulence model must obey. First 
of all, it is obvious that 

K> 0, £>0. (17) 


It can be shown that the conditions of realizability (17) are satisfied for homogeneous tur- 
bulent flows if: 


dl< ^ n d£ 

it- 0 


when 


K,£ 


0. 


Secondly, the turbulence model must be invariant under the Galilean transformation: 


x* = x — Uq£ 


(18) 


where Uo is any constant velocity. In this sense, the turbulence model must behave like 
the original Navier-Stokes equations which are Galilean invariant. Any model violating 
invariance under (18) is physically incorrect. 

Using the simplest closure for the deviatoric part of r,j: 

~K 2 

Tij — 2urSij = 2C )L — S{j (19) 

(where ~ 0.085 is a constant) (16) is transformed into a familiar form. It can be shown 
that (16) with (19) satisfies the realizability conditions 7T > 0 and £ > 0 provided that 
C e2 > 1 and Cfi > 0. These and other interesting properties of equations (16) and (19) have 
been recently reviewed by Speziale 8 . 

Let us introduce the dimensionless variables 


V = 


SK_ 

£ 5 


k- = JL 

I<o 


where S — (2 SijSij) 1 ^ 2 and Kq is the initial turbulent kinetic energy. In homogeneous shear 
flow where V 2 /lT = V 2 £ = 0, equations (16) and (19) have the simple form: 


dK‘ 

dt* 


= IC(C - V -') 


~ - l)r) 2 + (C £2 - 1 ) 


( 20 ) 
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where t* = St. Equation (20) has a single fixed point given by 


Vo 


Ce 2 — 1 


T 1/2 


C^Cex - 1) 


( 21 ) 


which is obtained by setting drj/dt = 0. This fixed point yields asymptotic solutions for K 
and £ of the form 

K cxe M \ £<xe Xt ' 
where the growth rate A is given by 8 

,_r C ll (C e 2 ~ C ei ) 2 I* 

\{C e \ - l)(C e2 - 1). 



Physical and numerical experiments on homogeneous shear flow - for not too large values of 
the dimensionless shear rate rj - indicate that indeed, rj — ► rjo and K,£ oc e Xt ’ for t* 1. 
This means that the solutions of any turbulence model must be attracted to this fixed point 
when rj is not too large. The fact that the simple K — £ model is capable of describing such 
a non-trivial behavior is remarkable. 

Another important consequence of equations (16) and (19) is the Reynolds number inde- 
pendence of the von Karman constant in turbulent channel flow. It is easy to show that the 
normalized dimensionless velocity profile U + is given in the region of constant energy K + 
by: 

U + = -In y + + B (23) 

K 

where k is the von Karman constant. The dimensionless variables y + and U + will be discussed 
in more detail in Section IV. It follows from the model (16) and (19) that the von Karman 
constant is independent of Re. This is the result of the cancellation of the singular terms 
assumed in the derivation of (13). If, on the other hand T\ + T 2 = 0(Re t ) and the 
assumption (13) is incorrect, then the von Karman constant k and the constant B in (23) 
will depend on Re t . This would have a significant impact on engineering calculations since 
(23) enters in the expression for the friction coefficient. Surprisingly, even today, there is 
no compelling experimental evidence for the constancy of B and k, although the prevailing 
opinion is that B and k are constants. 

Now let us discuss the closure assumption (19), which can be rewritten as: 


Tij = -2 CpKtjij (24) 

where rjij = SijK j£ is a dimensionless tensor which is equal to zero in isotropic turbulence. 
When rjij is small, the deviation from the isotropic solution is small. Unfortunately, in many 


5 


practical applications of the K — £ model, like channel flow, the values of parameter 77^ vary 
from T] ~ 20 in the buffer region to 77 ~ 3 in the logarithmic layer (where 77 = ||?7ij||). 

The dimensionless tensor 77 , j can be used for a perturbation expansion of r,y in powers 
of 77 when the departures from isotropy are small: 

T ij = -K[2C liVij + 0(r 1 l) + ...}. (25) 

In this case, the theory of strained turbulence can be formulated in terms of a double expan- 
sion in powers of the two dimensionless small parameters 7?e -1 and 77, j. This expansion gives 
rise to the anisotropic eddy viscosity explored by Speziale 9 and Rubenstein and Barton 10 . 
The existence of the second dimensionless expansion parameter 77 is not reflected in the 
functional form of the standard K — £ model. This may account for the relatively poor 
performance of the standard model in flows with large values of 77, j. 

It is apparent that a systematic derivation of turbulence models will enable us to address 
some of following questions: 

(i) Is it true that T\ + T2 = 0(1) as in (13)? 

(ii) Does the model (16) include all of the relevant effects or is something missing? 

(iii) Is it possible to derive a more complete Galilean invariant model than (16) satisfying 
the constraints discussed above, namely: 

(a) K > 0; £ > 0 (realizability) 

(b) dt]/dt = 0 when 77 — » 770 (fixed point) 

(c) K = 0(1)? 

Condition 3(a) and the requirement of Galilean invariance are the basic constraints, while 
constraints 3(b)-(c) are based on experimental data and must be accepted with some measure 
of skepticism until theoretically justified. 

II. Renormalization Group Methods and Turbulence Models 

The equations of motion (8)-(ll) describe general turbulent flow. We assume that the 
turbulence is driven by mean strains and decompose the velocity field into mean and fluc- 
tuating parts (v = U + u) as in Section I. We also make the basic assumption that the 
turbulence statistics are homogeneous when the mean strain Sij is not too large. Next, we 
assume that the fluctuating velocity field u is governed by the Navier-Stokes equations driven 
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by a random force, chosen in such a way that the global properties of the resulting field are 
the same as those in the flow driven by the mean strain Sij. 

To derive the equations for u, we consider an infinite domain in which a Newtonian fluid 
is stirred by a Gaussian random force. The force is defined by its correlation function in 
wavevector and frequency space (see Yakhot and Orszag 2 ), 




2Do(2Tr) d+l k y P(j(k)6(k + k')6(u> + u/), Al < k < oo 

(26) 


= 0, A s < k < A l 

where k is the wavevector, u is the frequency, k = |k| and d is the number of space di- 
mensions. The parameter y is chosen to give the Kolmogorov form of the small scale en- 
ergy spectrum; and in three dimensions y = d = 3. Statistical homogeneity in space and 
time is guaranteed through the factors £(k + k')(5(cj + u/) in (26). The projection operator 
Pij(k) — 6^ — kikj/k 2 makes the force statistically isotropic and divergence-free. In the limit 
of infinite Reynolds numbers, Al — ► 0 (but A l » A a — » 0). The forcing function (26) re- 
flects the fact that turbulence is usually driven by hydrodynamic instabilities with the most 
unstable mode at k — Al. The traditional view is that this primary instability generates 
velocity fluctuations with k > Al such that there is a direct energy cascade to small scales. 
The energy cascade may itself generate a universal random force of the form (26), but this 
has not yet been demonstrated rigorously. 

Our model is the forced Navier-Stokes equations: 


OV 1 

+ VjVjVi = -Vip + v 0 W 2 Vi + fi 


ViVi = 0 


(27) 


where the Gaussian force f is given by (26) and the density p has been absorbed in the 
pressure p. In contrast to the earlier development of RNG models in Ref. (2), the only new 
feature of the model (26)-(27) is the infrared cutoff of the random force: < fifj >= 0 when 
0 < k < Al- This property, which is usually unimportant, is needed if we wish to derive an 
equation for the mean rate of energy dissipation £. 

The dynamical equations for the total kinetic energy per unit mass K = and the 
homogeneous part of the instantaneous rate of energy dissipation per unit mass £ = u 0 ( V^Vj) 2 
are obtained from (27), 


dK 

dt 


+ ViViK = -£ + v 0 V 2 K - A7i{vip) + n,/,- 


d£ 

dt 


+ u,V t £ = 2^(Y J t> 1 )(Y J / t ) - 2t'o(V J u,)(V J t>,)(V / u,) 


-2i/$(V jV ivi) 2 - 2i/ 0 (VjU,)(V,Vjp) + i/ 0 V 2 £. 


(28) 

(29) 
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In the equation for £, driven by large-scale features such as boundary and initial con- 
ditions, T\ balances T \ to leading order. As we shall see, for steady-state flows (27)-(29) is 
sustained by the force f; T 2 is balanced by both T\ and the random-force contribution to the 
^-production given by Ps = 2z/ 0 (VjVi)(Vj/ t ). 

We seek equations for the mean values U =< v >, K =< K > and £ =< £ >, averaged 
over an ensemble of the fluctuating velocity field. To find these equations, we shall use the 
dynamic renormalization group and the ^-expansion. Since the renormalized equations for 
K and £ may not be trivially related to the renormalized Navier-Stokes equations, the RNG 
procedure must be applied to all of the equations (27)-(29). 

The renormalization group applied to (28)-(29) is described in detail in Yakhot and 
Smith 11 . Here we present the main results. The equations of motion averaged over the 
fluctuating velocity field v are: 


dUj 

dt 


+ U • W t = -V iP + 





d£ 

dt 


dK TT 0 d ( dl<\ 

dt + U ‘ VA ” s TijSij + dxi dxi ) 

+ u • v? = -Ce^ iiSii - - n + A 


where 


71 — 2 is 0 Si 


~5ui~5ui 

dxi dxj ’ 


(30) 


(31) 


Cc x — 1.42, Ce 2 = 1.68, and ax = ot£ = 1.39. The relations (30)-(31) are obtained in the 
lowest order of the effective Reynolds number: 



2£> 0 

i/A? 


oc e ~ 0(1) 


(32) 


where A , ~ 2i t/ L and v is the effective viscosity (7). The parameter e = 7 — y is characteristic 
of the random noise driving this artificial flow. It was also shown 11 that in this order of 
the expansion, T\ + T 2 = 0(1) which confirms the hypothesis leading to (13). In sheared 
turbulence where Sij / 0, the Reynolds decomposition generates contributions of the order 


c dui dui 
U ° ij dxi dxi 


(33) 


The isotropic and homogeneous random velocity field v was assumed to be governed 
by the equations (2)- (4). The renormalization group and the e-expansion were applied 
to evaluate (33). This systemmatic procedure generated an expansion in powers of the 
dimensionless parameter rjij = SijK/S. The contribution to 7Z in (31) is small for weakly 
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strained turbulence and large in the rapid distortion limit when 77 — > 00 . All attempts to 
close (31) using the methods based on the e-expansion were fruitless. In the next section we 
shall show how to evaluate (31) by applying the general criteria outlined in Section I. 


III. Generalized K — £ model 


Neglecting 7Z in (31) and using the low order closure (19) for r,j, we obtain equations of 
motion, which are those of the standard K — £ model: 


DK , 7 d 

— = 2vrSij - £ + ^ 


OIK V 


dK ' 

dx ij 


(34) 


— = 2 C e , w u T S, 


u 


C £ 2= + £- 

I< axi 


Ot£ v 


d£T 

dxi , 


where D/Dt = d/dt + U • V. 

This model satisfies all of the criteria formulated in Section I. In homogeneous shear flow 

~ 4.38 and 


Cfo- 1 


where dUi/dxj = S6n8j2, equations (34) have the fixed point T] 0 — [c^cd- 1 ) 
for St > 1 the kinetic energy and dissipation grow exponentially: 


1/2 


K £ Af . 

=- oc — oc e M 
I< 0 £ 0 

with a growth rate A = 0.142, which is very close to that obtained from numerical and 
physical experiments. These results will be discussed in more detail in Section V . 

The model ( 34 ) is based on the assumption that 77 is small; however, the numerical value 
of 77 4.38 does not satisfy this constraint. Furthermore, in the logarithmic region of a 

turbulent channel flow, the value of the dimensionless shear 77 = SK j £ — 3 is also not small. 
This means that the contribution to Eq. (30) due to 7 Z cannot be neglected there. 

Iterating the expression for 7 Z using the Navier-Stokes equations will generate a power 


series: 


' SI<' 


(35) 


7Z — vtS 3 ^2 r n , T 

71=0 \ c / 

where S = (‘2S lJ S lJ ) 1 ^. It is clear that when S — > 0, 1Z goes to zero faster than j^vjS 1 and 
may be neglected in (30). To use (35) in the case of large shear, let us consider as a model 
only the subset of terms corresponding to the geometric series 


= v T s 3 £(-/?)’ 


’SK 


3 n 


n=0 


vtS 3 

1 + fir) 3 * 


(36) 
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We do not know the values of the coefficients r n and cannot evaluate the magnitude of 
71 — 71°. However, this form of 7Z° is plausable since (36) implies that 


71° = 0 


for all values of rj. Thus, 71° vanishes when £ — * 0, a constraint that also follows from 
the definition of 7Z given in (31), so that (36) maintains the satisfaction of the realizability 
conditions (17). Assuming that the fixed point value T] o ~ 4.38, obtained in the limit of 
small Sij, is invariant to dropping all terms but those in (36) we postulate that: 



71 = 


»^ 3 (i - j) 


i + /V 

Equation (37) can be rearranged into the alternative form 


(37) 


71 = 


(38) 


1 + /V I< 

from which it is clear that the crucial constraint that 71 —* 0 as £ — * 0 is satisfied. The 
limits of weak and strong strains can also be easily discerned from (38). It is clear that in 
the limit of weak strains where rj — > 0: 


r 


n ~ o. 


oo: 


On the other hand, in the limit of strong strains where rj 

P Vo I< 

which becomes singular. From the definition of 71 given in (31), it follows that 


(39) 


(40) 


£ 2 

ft = 2 vA l} s;-= 


(41) 


where 


dui dui 
A 0 dxi dxj 

Aij = - 


Q* — lJ 
Zij- s • 


oo. 


It is a simple matter to show that ||Aij|| 2 < | and ||5 *j|| 2 < \ so that in the limit as rj 
it follows that 7Z = 0(v) consistent with (40). Furthermore, in the limit as £ — * 0, 71 — ♦ 0 
as 0(£ 2 ) and in the limit as v ~ * 0, 71 — > 0 as 0(v 2 ) or 0( v 3 ) (which follows from local 
isotropy; c.f. Monin and Yaglom 12 ). These results are generally consistent with (39). 

It is clear that the model (38), while not rigorously derived by RNG methods, satisfies 
a variety of important consistency conditions that place severe constraints on the possible 
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functional form of K. In particular, no finite number of terms in the series (35) can give the 
desired asymptotics. While the analyticity of 7 ^( 77 ) in the limit as rj — > oo is not necessary 
we choose to assume it for simplicity. If this is the case, then the postulate (38) for 7Z is 
quite reasonable. 


IV. The von Karman Constant 


The proposed model for TZ- given in (38) contains one undetermined constant /? which will 
now be related to the von Karman constant. In order to do this, we rewrite the dissipation 
rate transport equation in the equivalent form 


ft + u ■ V£ - - c *h iSii ~ + h ( ae,, 3 (42) 


where the coefficient C| 2 is given by 


C * 2 = 1.68 + 


O 3 0 - £) 

l + ,v 


(43) 


In the log layer of a turbulent boundary layer the turbulence production T{j $ij — £ and, 
hence, 77 = Ify/C^ which renders C7 | 2 constant. Furthermore, convection effects (i.e., DS/Dt ) 
can be neglected and the mean velocity assumes the logarithmic form (23) where 


y 


+ _ 


yUr 

V 



U T 


and u T is the shear velocity. With these assumptions, it is a straightforward calculation to 
show that (42) yields the constraint 13 


K — 


{Pet ~~ 


as 


(44) 


which establishes a direct relationship between the von Karman constant k and the constant 
l 3 in the model (38) for 7 Z. It follows that when 


0 = 0.012 

the von Karman constant k = 0.4. Fortunately, k is not very sensitive to variations in 0 (in 
applications, 0 can be varied between 0.01 and 0.015 with only small changes in n). With 
this selection of /3 (and the fixed point T} 0 = 4.38) we now have a closed model for U. Since 
the model yields a log layer with a von Karman constant k of 0.4, it will give good results 
for turbulent channel flow. The model will be subjected to a more severe wall-bounded flow 
test in Section VI where turbulent flow over a backward facing step will be considered. 
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V. Homogeneous Shear Flow and the Relaxation Time Approximation 


In this section we report tests of the new model in homogeneous shear flow where an 
initially isotropic turbulence (with turbulent kinetic energy K 0 and dissipation rate £ 0 ) is, 
at time t = 0, subjected to a constant shear rate S with the corresponding mean velocity 
gradient tensor 


dUi 

dxj 


— S8i\6j2- 


The model given by equations (19), (30), and (38) yields a simple set of coupled nonlinear 
ordinary differential equations since I< = K(t) and £ — H(t). The solution obtained for the 
time evolution of the turbulent kinetic energy is compared in Figure 1 with the large-eddy 
simulation of Bardina et al. 14 Despite the fact that the model yields an excellent long-time 
growth rate for K (i.e., the model predicts that A = 0.142 in comparison to physical and 
numerical experiments 15-16 which indicate that A is in the range of 0.14 - 0.16), the early time 
values of K are overpredicted. This results from the use of an eddy viscosity representation 


for T{j which responds too quickly to the application of the shear at t = 0. 

The early time behavior can be described more accurately when a relaxation time ap- 
proximation is introduced. In this approximation, we allow the Reynolds stress r,j to relax 
to the eddy viscosity model (19) as follows: 


— (r.j + 2 1'jSij) = -C TR j=(Tij + 2 v T Sij) (45) 

for homogeneous flows where Cj R is a relaxation coefficient. It is a simple matter to show 
that as a direct consequence of (45) 

it = ~ Ct »T T< i ~ 2C - {Ct - + c ‘ 2 - 2)Ys ‘i + («) 

When an initially isotropic turbulence is strained, Crow 17 showed that 

C„(C T „ + C„- 2) = ^ 

which yields Cj R ~ 3.5. Eq. (46) gives a simple Reynolds stress transport model. 

In addition to (45) for the relaxation time approximation, vtS 2 should be replaced by 
the production P = —TijSij in Eq. (37) giving the relaxation model 


n = 


Mi z£_Z 

1 + /V K 


(47) 


where for homogeneous shear flow P = — T 12 S. When the eddy viscosity model (in which 
P — 2vTSijSij) is introduced in (47), Eq. (37) is recovered. 
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In Figure 2, the time evolution of the turbulent kinetic energy in homogeneous shear 
flow obtained from the relaxation model (46) - (47) is shown. These results are in excellent 
agreement with the large-eddy simulation of Bardina et al. 14 This is quite encouraging con- 
sidering the fact that this relaxation model was calibrated independent of homogenous shear 
flow (i.e., it was calibrated by setting the value of the von Karman constant to 0.4 and by 
using the Crow result for plane strain turbulence). 


VI. Turbulent Flow over a Backward Facing Step 


Turbulent flow over a backward facing step is one of the standard test cases used to 
evaluate the performance of turbulence models in separated flows. The Kim, Kline and 
Johnston 18 experimental test case will be considered here for the backstep problem. For this 
flow configuration, the expansion ratio (step height: outlet channel height) E is 1:3 and the 
Reynolds number Re = 132,000 based on the inlet centerline mean velocity and the outlet 
channel height (see Figure 3). We will use the data for this test case as updated by Eaton 
and Johnston 19 for the 1980-81 AFOSR-HTTM Stanford Conference on Turbulence. Since 
the turbulence is statistically steady, we can use the simpler form of the model given by (37) 
(i.e., there is no need to incorporate relaxational effects). The mean turbulence equations 
are solved numerically by a finite volume method on a 200 x 100 variable mesh (see Thangam 
and Speziale 20 for more details on the numerical method). Inlet profiles for U, K, and £ are 
specified five step heights upstream of the step corner (a developing turbulent channel flow 
calculation is carried out that matches the experimental inlet mean velocity profiles of Eaton 
and Johnston 19 ); conservative extrapolated outflow conditions are applied thirty step heights 
downstream of the step corner. Along the upper and lower walls, the law of the wall (23) with 
k — 0.4 and B — 5 is used along with standard log layer boundary conditions for the kinetic 
energy and dissipation obtained from a production-equals-dissipation equilibrium hypothesis 
and the constraint of vanishing normal derivative of K at the wall. It was demonstrated 
recently by Thangam and Speziale 20 that, provided the flow field is well resolved, the use of 
these law of the wall boundary conditions leads to only minimal errors. 

Two sets of results will be presented based on the algebraic model (37) for 7Z: one set 
obtained by using the isotropic eddy viscosity model (19) and another obtained from the 
anisotropic eddy viscosity model of Speziale 9 given by 
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where Sij is the Oldroyd derivative of S tJ with convective effects neglected; C'd is a dimen- 
sionless constant which assumes a value of 1.68. For thin shear layers, this model is very 
similar to the anisotropic eddy viscosity model derived by Rubenstein and Barton 10 based 
on RNG methods. This type of anisotropic eddy viscosity model is obtained when terms of 
order 77? are maintained in the expansion (25) for r,j. 

In Figures 4(a)-(b), the computed mean velocity streamlines and mean velocity profiles 
obtained using the isotropic eddy viscosity model (19) are compared with the experimental 
data. 19 Reattachment is predicted at xr/ H ~ 6.6, a result that is within 6% of the ex- 
perimentally measured reattachment point xr/H 7.0. The computed turbulence intensity 
and turbulence shear stress profiles are compared with experimental data in Figures 5(a)-(b). 
The agreement is comparably good. 

By combining the model (37) for 1Z with the anisotropic eddy viscosity model (48), even 
better results are obtained. In Figures 6(a)-(b), the predicted mean velocity streamlines 
and mean velocity profiles are shown to compare extremely favorably with the experimental 
data. 19 The model predicts reattachment at xr/H 7.0 which is essentially the same as 
the experimental result. Furthermore, the model predicts a noticeable secondary separation 
bubble below the corner of the step consistent with experimental observations for this back- 
step flow. The agreement between the model predictions and the experimental data for the 
turbulence intensity and turbulence shear stress profiles is also excellent as shown in Figures 
7(a)-(b). The quality of these predictions is quite remarkable for a two-equation model. 

VII. Concluding Remarks 

The renormalization group formalism of Yakhot & Orszag for the development of turbu- 
lence models has been supplemented with scale expansions for the Reynolds stress and pro- 
duction of dissipation terms. Here, the extra expansion parameter is taken to be rj = SK/€ 
which is ratio of the turbulent to mean strain time scale. For the Reynolds stress, this 
approach - which leads to the development of anisotropic eddy viscosity models as well 
as Reynolds stress transport models - is analogous to that introduced by Rubenstein and 
Barton. 10 However, the present method is completely new for the modelling of the production 
of dissipation term 1Z which is neglected in most of the commonly used turbulence models. 
The interesting result for 1Z is that no finite order truncation of the expansion satisfies the 
necessary physical constraints; terms of all orders in the expansion must be retained. This 
complication eliminates the possibility of determining the model explicitly in closed form. 
However, a highly plausible form, with only one undetermined constant, is postulated here 
which satisfies all of the necessary physical constraints (i.e., realizability and consistency with 
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the weak and strong strain limits). The constant is calibrated by setting the von Karman 
constant k to 0.4. 

The new models have been tested for homogeneous shear flow and for flow over a back- 
ward facing step. Excellent results are obtained in both cases. For the case of homogeneous 
shear flow, the best results are obtained from the Reynolds stress transport model (i.e., the 
relaxational model discussed in Section V). On the other hand, excellent results are ob- 
tained with eddy viscosity models for the backstep problem. In all of these calculations, no 
ad hoc adjustments of the constants are made. The applications considered in the paper are 
restricted to simple shear flows since the current version of the dissipation rate transport 
equation has only been modelled to account for the effects of irrotational strains. Incorpora- 
tion of the effects of rotational strains, which can be important in turbulent flows involving 
curvature or a system rotation, is a difficult task that has not yet been achieved. The reduc- 
tion in the energy cascade that occurs in rotating isotropic turbulence and the stabilizing 
or destabilizing effects of rotations on homogeneous shear flow are but two examples. 
This more difficult problem of accounting for rotational strains using a comparable scale 
expansion technique will be the subject of a future study. 
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Time evolution of the turbulent kinetic energy In homogeneous shear flow. 
Relaxation model; o Large-eddy simulation of Bardina etal for e 0 /SK 0 - 0.296 


Figure 2 
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Turbulent flow over a backward facing step 
Figure 3 
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BACKWARD-FACING STEP: 




(b) Dimensionless mean velocity profile 

( Computations with isotropic eddy viscosity; 

o Experiments of Kim et al, 1980; Eaton & Johnston, 1981) 


Computed mean flowfield for the new RNG K-e model 
[E = 1 :3; Re = 132,000; 200x100 mesh] 

Figure 4 
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BACKWARD-FACING STEP: 



(a) Turbulence Intensity 
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(b) Turbulence shear stress 

Computed turbulence stressesfor the new RNG AC-emodel [E = 1 :3; Re = 1 32,000; 

200x100 mesh; computations with isotropic eddy viscosity; 

o experiments of Kim et al, 1980; Eaton & Johnston, 1981] 

Figure 5 
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BACKWARD-FACING STEP: 



X/H 

(a) Streamlines 
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(b) Dimensionless mean velocity profile 

( Computations with anisotropic eddy viscosity; 

o Experiments of Kim et a/, 1980; Eaton & Johnston, 1981) 


Computed mean flowfield for the new RNG fC-e model 
[E = 1:3; Re = 132,000; 200x100 mesh] 

Figure 6 
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BACKWARD-FACING STEP: 
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(a) Turbulence Intensity 
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(b) Turbulence shear stress 

ComputedturbulencestressesforthenewRNGK-emodel[E=1:3;Re = 132,000; 

200x100 mesh; computations with anisotropic eddy viscosity; 

o experiments of Kim et al, 1980; Eaton & Johnston, 1981] 

Figure 7 
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